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We study how parameter error forecasts for tomographic cosmic shear observations are affected 
by sky coverage, density of source galaxies, inclusion of CMB experiments, simultaneou fitting of 
non-dark energy parameters, and the parametrization of the history of the dark energy equation- 
of-state parameter w(z). We find tomographic shear-shear power spectra on large angular scales 
(I < 1000) inferred from all-sky observations, in combination with Planck, can achieve u(wo) = 0.06 
and o(w a ) = 0.09 assuming the equation-of-state parameter is given by w(z) = wo + w a (l — a(z)) 
and that nine other matter content and primordial power spectrum parameters are simultaneously 
fit. Taking parameters other than wq, w a and fi m to be completely fixed by the CMB we find errors 
on wq and w a that are only 10% and 30% better respectively, justifying this common simplifying 
assumption. We also study 'dark energy tomography': reconstuction of w(z) assumed to be constant 
within each of five independent redshift bins. With smaller-scale information included by use of the 
Jain & Taylor ratio statistic we find o~{wi) < 0.1 for all five redshift bins and o~(wi) < 0.02 for 
both bins at z < 0.8. Finally, addition of cosmic shear can also reduce errors on quantities already 
determined well by the CMB. We find the sum of neutrino masses can be determined to ±0.013 
eV and that the primordial power specrum power-law index, ns, as well as dns /dink, can be 
determined more than a factor of two better than by Planck alone. These improvements may be 
highly valuable since the lower bound on the sum of neutrino masses is 0.06 eV as inferred from 
atmospheric neutrino oscillations, and slow-roll models of inflation predict non-zero dns /dink at 
the forecasted error levels when \ns — 1| > 0.04. 

PACS numbers: 98.65.Dx,98.70.Vc,95.35.+d,98.80.Es 



I. INTRODUCTION 

The great mystery of twenty-first century cosmology 
is the dark energy — a smooth, negative-pressure com- 
ponent evoked to bring the Universe to critical density 
( e -g-iBI3) an d t° dim high-redshift supernovae 0,0, Of. 
Observers are constraining, or planning to constrain, this 
curious component in a number of ways, both through its 
influence on structure formation and on geometry. Here 
we consider how well the component can be probed by to- 
mographic cosmic shear observations of large-scale struc- 
ture. 

Cosmic shear is a very promising technique for study- 
ing the large-scale structure and constraining the cos- 
mological parameters that govern its evolution. For a 
recent review see @- The shear signal only depends on 
the mass distribution; it is immune to uncertainties in 
the relationship between mass and light that plague other 
observations. On large scales, given a model, the statisti- 
cal properties of this mass distribution can be calculated 
with very high accuracy, thus there is a well-understood 
relation between theory and observable. 

Cosmic shear also naturally probes the chief redshift 
range of interest for constraining dark energy: < z < 2. 
In this redshift range, dark energy becomes the dominant 
contributor to the energy density and thus has significant 
effect on the expansion rate, H(z). The expansion rate 
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affects the shear power spectrum both through its effect 
on the growth rate and also by altering how length scales 
at a given redshift project into angular scales. At z > 2, 
for the simplest models, the dark energy is highly sub- 
dominant and has little influence on structure formation. 

Multi-band photometry makes possible tomographic 
cosmic shear, in which source galaxies are binned accord- 
ing to their photometrically-determined redshifts, and 
separate shear maps are made for each source galaxy red- 
shift bin. This technique allows one to study the evolu- 
tion of the density field as a function of redshift, restor- 
ing the information otherwise lost by projection over the 
broad lensing kernel. 

Many papers 0, H, S E3, E3, E3, E3, II^I have stud- 
ied how well tomographic shear correlations can be used 
to constrain dark energy parameters. All but one of these 
error forecasts were done in the limit of other cosmo- 
logical parameters being perfectly known a priori. This 
assumption is justified if CMB (or other) observations 
can be used to determine the parameters with negligibly 
small errors. We investigate this assumption by letting 
eight other cosmological parameters float, constrained 
only by the cosmic shear data themselves and various 
CMB observations. For WMAP, Planck, and a future 
polarization satellite described below we find the con- 
straints on the dark energy parameters degrade by 70%, 
30% and 4% respectively. 

We parametrize the time-dependence of w in two differ- 
ent ways. First, following Linder |16| . we write w(z) — 
wq + w a (l — a(z)). This parametrization has the ad- 
vantage of tending toward a constant at high redshift, 
preventing high-redshift dark energy dominance that can 
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occur with the parameterization w{z) = wq + w z z. It is 
therefore a convenient parametrization when combining 
z < 2 observations with z ~ 1100 CMB observations. 

Weller & Albrecht have argued that all mod- 
els in the literature can have their luminosity distances 
as a function of redshift fit by the w(z) = wq + w z z 
parametrization to 1% accuracy in the < z < 2 range. 
However, dark energy is so poorly understood that it is 
worth considering w(z) histories that do not occur in any 
existing scalar field quintessence models. Therefore we 
also study a less model-dependent parameterization, one 
in which w(z) is constant within each of five independent 
redshift bins. An approach somewhat similar to ours 
is that of |l8fl who parametrized w(z) with its principal 
components, given supernovae data as expected from the 
Supernova/ Acceleration Probe (SNAP). We have found 
that data that give similar constraints in the wq , w a plane 
may give strikingly different constraints in the space of 
redshift bins of w. 

In this paper we concentrate on what is possible with 
the data on large angular scales where theoretical inter- 
pretation is most straightforward. There is also great 
potential for data on smaller angular scales although the 
state of theory must be improved to take full advantage 
of it. Potentially great use can be made of two and 
three-point functions out to I = 10 5 (e.g., 0). Even 
at I < 3000, including the 3-point function can improve 
errors by a factor of 3 ■ Counting of mass clusters is 
another promising way to exploit the small angular scale 
data 0. 

We do consider one method of usin g th e small angular 
scale data here. Recently [2(J and |2l| have proposed 
using the data in such a way as to be insensitive to 
the uncertain fluctuation statistics and only dependent 
on easily-calculated geometric quantities. We combine 
a conservative use of the ratio statistic of [2(J with the 
large angular scale shear-shear power spectra. We find 
that the addition of the ratio statistic greatly improves 
our reconstruction of w(z) even with a conservative fore- 
cast of how well the ratio statistic can be measured. 

In addition to dark energy, another exciting applica- 
tion of high-precision cosmic shear data is to determine 
the mass of the relic neutrinos. In the first cosmic shear 
+ CMB forecasting papers, jU it was shown that 
combining Planck and all-sky cosmic shear could deter- 
mine the sum of the masses of the different flavors of 
relic neutrinos to 0.02 eV with w fixed to -1. Recently 
[24| showed that a post-Planck all-sky CMB polariza- 
tion mission could be used to determine the sum to 0.04 
eV without cosmic shear data. The sum is at least 0.06 
eV in order to explain atmospheric neutrino oscillations 
[2q hence a detection is assured, though possibly one 
of weak statistical significance. Here we show that all- 
sky CMB polarization observations combined with all- 
sky weak lensing observations can reduce this error to 
0.027 eV or, with the inclusion of the ratio statistic, to 
0.013 eV even with simultaneous determination of wq and 
w a . This reduction in error may change a two a detection 



to a five a detection. 

We discuss the shear-shear correlations in section II 
and the binning of the dark energy history in section III. 
In section IV we discuss the parameters of the future 
weak lensing and CMB experiments we consider and in 
section V outline our error-forecasting method. In sec- 
tion VI we give results for constraints on the dark energy 
and neutrino mass and all the other cosmological param- 
eters we assume our simultaneously fit. In section VII 
we see how the results improve with the addition of the 
ratio statistic. 



II. TOMOGRAPHIC COSMIC SHEAR-SHEAR 
CORRELATIONS 

The size and shape of galaxies is altered by gravita- 
tional lensing. The effect of lensing is described by the 
2x2 transformation matrix Aij given by 

_ d9\ _ ( 1 - k - 7i 72 \ m 
y -d&~\ 72 1 - k + 7i J ■ W 

where 9 S is the angular location in the source plane for a 
light ray appearing at 9 in the image plane. The conver- 
gence k describes the magnification or demagnification 
while the shear components, 71 and 72, quantify the dis- 
tortion of the shape |Ml27| . 

The shear components can be inferred from measure- 
ment of galaxy ellipticities. Galaxies have their own 
intrinsic ellipticity and an additional ellipticity due to 
lensng so that 

e" = e^ int + e"i en8 (2) 

where /i labels the galaxy. The average shear in a pixel on 
the sky can be estimated from averaging over the galaxies 
in the pixel: 

ti = 1 Vc' 1 = 1 V (e^ + ) 

' AT ■ ' ' + OAT ■ ' ' V +.int +.lons/ 

^ iv pix ^ Jv pix 

^ 2 = 9~/V 5Z C x = OAT . , int + e x , Ions) (^) 

fj, n 

where + and x denote the two orthogonal ellipicity com- 
ponents. In the absence of correlations between the in- 
trinsic ellipticities, the rms error in the measurement of 
each shear component is 

O (71 ) = a (72 ) = 7rms / \/N^ (4) 

where j Tms is the rms intrinsic shear of the galaxies and 
Npi x the number of galaxies in the pixel. 

Maps of the shear components can be decomposed into 
even parity E modes and odd parity B modes, just as is 
the case with linear polarization Stokes parameters, Q 
and U psl |29| . Assuming uniform coverage, the noise in 
the E modes will have a diagonal covariance matrix N 
given by 

Nl m B,l'm'B' = N] 3 8w 8 mm '&BB' ■ (5) 
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FIG. 1: The shear-shear auto power spectra. The 8 solid 
curves are the shear power spectra from each of the galaxy 
source planes, Bi to B$. Dotted curves are the linear pertur- 
bation theory approximation. From bottom to top the source 
plane redshift ranges are By. z G [0.0,0.4], B 2 : z G [0.4,0.8], 
B 3 : z G [0.8,1.2], B 4 : z 6 [1.2,1.6], B 5 : z G [1.6,2.0], B e : 
z G [2.0,2.4], B 7 \ z G [2.4,2.8] and B s : z G [2.8,3.2]. The er- 
ror boxes are forecasts for G4ir (see Table I). The top dashed 
curve is the shear power spectrum for the CMB source plane. 
The error boxes are forecasts for CMBpol (see Table II). 



where 



N? = jLMN^ ( fi ) 



and B labels the source redshift bin. The signal con- 
tribution to the covariance of the shear E modes, S, is 
correlated across bins, but otherwise diagonal: 



Sl m B,l'm'B' — C BB Sw5 mm '. 

The shear angular power spectra are given by 

T 2; 
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BB' 
I 



= ir / drrwB ( r ) wB '( r ) A l( k >> 



(8) 



where the window function for sources in redshift bin, B, 
is 



We use Peacock and Dodds' formulation [30J to cal- 
culate the non-linear density power spectrum, A|(fc,?*). 
We only use data with I < 1000 to reduce our sensitiv- 
ity to errors in the non-linear corrections and departures 
from Gaussianity |3lj| . 

We assume that photometrically-determined redshifts 
can be used to sort the source galaxies into eight redshift 
bins with Az = 0.4 from z = 0.0 to z — 3.2. The shear 
power spectra for these eight source bins are shown in 
Fig. 1. The higher the source bin redshift, the higher the 
shear power. The shear power spectra are related to the 
deflection angle power spectra by 



1(1 + l)Cr/(2n) 



(11) 



thus the area under the curves in Fig. 1 gives the contri- 
bution to the deflection angle variance from a given range 
of I. The expected errors on the auto power spectra are 
given by 



AC BB = 



(21 + l)f sky 



BB 



N, 



BB \ 



(12) 



Shear maps can also be inferred from observations 
of CMB temperature and polarization. The top curve 
in Fig. 1 is the shear power spectrum for sources at 
z ~ 1100. The auto and cross-correlation (with the eight 
other shear maps) power spectra for the CMB shear map 
is also given by Eq. |H1 However, the window function, 
due to the small thickness of the last-scattering surface 
is simply given by 



W 



CMB 



r(z)] 



r(z s ) - r(z) 
r(z s )r(z) 



(13) 



where z s is the redshift of the last-scattering surface, de- 
fined as the peak of the visibility function. We calcu- 
late the error in the reconstructed shear maps using the 
quadratic estimator of 32]. The noise level in Fig. 1 
is for the CMBpol experiment described in section IV. 
At these noise levels, a maximum-likelihood shear map 
would have up to two times lower reconstruction noise 
than the quadratic estimator shear map |33| . 

Thus we have grouped the data into nine different red- 
shift bins for lensing tomography: eight bins with galax- 
ies as the sources (from z = 0.0 to z = 3.2) and a ninth 
bin on the CMB last-scattering surface. 



III. BINNING DARK ENERGY HISTORY 



1 2 r^max) (r 1 - r) 

W B (r) = / dr'± -N B (z')dz'/dr'e(r'-r), of state parameter as 

n B r J r ( z B ) r' 

and the average number density of galaxies in redshift 
bin B is 



We parametrize the history of the dark energy equation 
as 

(z)=J2Mz)w l (14) 



n B 



dr'N B (z')dz' /dr' 



where 
( 10 ) X, 



otherwise. 



(15) 
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We define 5 bins: 4 bins with equal redshift spacing Az — 
0.4 and the last bin from z = 1.6 to the last scattering 
surface. We choose this binning because the data are 
sufficiently powerful to place interesting constraints on 
this space. Finer binning, particularly beyond z = 1.6, 
would lead to very strong degeneracies between the bins. 

Note that the discrete parametrization of the dark en- 
ergy equation-of-state parameter is well-matched with 
what we observe in weak lensing experiments with binned 
source planes. The variation of lUj in bin % influences the 
dark energy density in all the lower redshift bins. Like- 
wise, the shear map from galaxies in bin i are affected 
by the large-scale structure in those same lower redshift 
bins. 

This discrete binning of w(z) does lead to complica- 
tions on very large scales where dark energy density fluc- 
tuations can be important. While it is straightforward 
to calculate the expansion rate and therefore distance as 
a function of redshift for w(z) given by Ea. (|14l) . it is not 
clear how to calculate the fluctuations in the dark energy 
density and pressure for a model with non-continuous 
w(z). On small scales, the dark energy fluctuations (at 
least for the simplest scalar-field models) are damped 
out. The critical scale Uq (Iq in harmonic space) is calcu- 
lable from the coupled differential equations of the matter 
and the dark energy fluctuations: 

5 m + 2HS m = ATrG(p m 5 m + Sp Q + i5p Q ) 7 (16) 

S Q + 3H5 Q + (k 2 + 5q = S m [(1 + w Q )p Q ] 1/2 

where S m = 6p m / Pm, Sq = Spq/pq and a dot denotes 
derivative with respect to proper time. In the limit where 
the inhomogeneity in Q-fields is negligible, the evolution 
of S m is decoupled from the dark energy fluctuations yj, 
l35| . In this limit, the influence of dark energy is entirely 
through H(z). 

At k < ^/\d 2 V/dQ 2 \, Sq grows via gravitational in- 
stability and will cluster like a matter component. The 
evolution of 5 m will be affected by Sq and will therefore 
be &;-dependcnt in the era of dark energy domination. At 
k > y/\d 2 V / dQ 2 \, SQ is damped out and the evolution 
of S m will be independent of SQ. The characteristic scale 
kQ of this regime is given by 

k Q = o.ooi^(i -„) (^^- an+ t; a ^ )' (17) 

where kQ is in the unit of /iMpc -1 [sfij . 

The characteristic angular wavenumber is given by 
Iq ~ kQr(rf) where r is the angular diameter distance 
and r\ is conformal time. Assuming our fiducial model, 
Iq is around 20 for the CMB shear map and even less 
for the galaxy source plane shear maps. We therefore 
restrict our analyses to I > 40 where the shear maps will 
be unaffected by dark energy fluctuations. 

Analysis of real data could proceed in two stages. At 
the first stage, one would fit the data with the modes I > 
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TABLE I: Weak lensing experimental parameters assumed. 
'G' and 'S' are for ground- and space-based observations re- 
spectively which have different source redshift distributions. 
Units for the total source sky density, n to t, are 1/arcmin 2 and 
the per-component rms intrinsic shear, 7 rms , is evaluated at 
z = 1. 
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TABLE II: CMB experimental parameters assumed. Planck 
has more channels which we assume are used entirely for con- 
trol of foregrounds. 



40 where the effect of SQ on matter density perturbations 
is negligible. With the best-fit Wi determined, one could 
use interpolation to generate a smooth w(z) as well as 
the first few eigenmodes of fluctuations about w{z). The 
amplitudes of these eigenmodes could then be fit to all 
the data. 



IV. EXPERIMENTS 

Weak lensing surveys are beginning to make their im- 
pact on cosmology 6j. For example, Contaldi et al.[37j 
find the combination of early Red-Sequence Cluster Sur- 
vey (RCS) data and WMAP data give cr 8 = 
0.89 ± 0.05 and n m = 0.30 ± 0.03. Ongoing surveys of 
tens to hundreds of square degrees include the Deep Lens 
Survey (DLSjlEl), the Canada- France-Hawaii Telescope 
Legacy Survey (CFHTLS jz^hthe National Optical As- 
tronomy Observatory (NOAOjl^) Deep Wide-Field Sur- 
vey, and the Red-sequence Cluster Survey (RCS jzH)- 
Longer term there are projects that may cover 1000 
square degrees or more. These are the SuperNova Accel- 
eration Probe (SNAP HjJ), Pan-STARRS H3 and the 
Large-aperture Synoptic Survey Telescope (LSST 82]). 

We consider two sets of reference surveys, ground- 
based and space-based which have different total source 
density and source redshift distributions. For the ground- 
based surveys we use a galaxy redshift distribution for a 
limiting magnitude in R of 26 inferred from observations 
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with the Subaru telescope |39j ■ The shape of this distri- 
bution is well-described by the following analytic form: 



dn/dz oc z exp — (z/1.2) 
dn/dz oc z' exp —(z/1.2) 



for z < 1 
for z > 1. (18) 



We use this distribution with the modification that half 
of the galaxies in the 1.2 < z < 2.5 range are discarded 
for insufficiently accurate photometric redshifts. The am- 
plitude of the distribution is such that, after this cut, the 
number density of galaxies is 65 per sq. arcmin|8^|. 

Space-based surveys have the advantage of no sky 
noise, allowing them to reach higher limiting magnitudes, 
and the possibility of additional bands in wavelength 
ranges for which the atmosphere is opaque. The latter 
advantage means photometric redshifts can be done even 
in the problematic 1.2 < z < 2.5 range 40]. For the 
space-based surveys we assume 



dn/dz 



' exp (—2/1.5) 



(19) 



which was also assumed by , 21] and which is a close ap- 
proximation to what was used in a study of the SNAP 
300 sq. degree survey |4jj ■ This study predicted a num- 
ber density of ^100 galaxies per sq. arcminutc. 

The disadvantage of space is that the telescopes are 
smaller due to the high cost of putting large mirrors in 
space. This means that ground-based surveys can cover 
much more sky. Our two ground-based surveys are de- 
noted by G4-7T (an all-sky survey) and G2-7T (a half-sky 
survey) and our two space-based surveys are called SI 000 
and S300 with the numbers denoting the observed area 
in square degrees. Their properties are summarized in 
Table I. Baseline plans for LSST include a G2-7T survey. 

An important quanity for forecasting parameter con- 
straints from any cosmic shear observation is 7rms . 
This quantity is potentially experiment-dependent, with 
ground-based observations having a higher effective 7 rms 
than space-based ones j4?J ^J. But for the low- noise 
ground-based observations we consider here, 7 rms may be 
just as goo d as from space. In fact, ground-based obser- 
vations ^2 have achieved very low 7 rms . Jarvis et al. 
find for the lowest-redshift galaxies in their 75 sq. degree 
CTIO survey that the shape noise contribution to the 
per-component rms shear, 7 rms , is 0.15. The shape noise, 
as inferred from these data, increases approximately lin- 
early to 0.22 by z = 20- This trend with redshift is also 
consistent with simulations of SNAP weak lensing obser- 
vations, based on HDF data, that show 7 rms , averaged 
over all galaxies suitable for cllipticity measurement, in- 
creasing from 0.21 to 0.25 with increasing observing time 
per unit solid angle [4(j. We therefore model 7 rms as 



7rms (z) = 0.15 + 0.0352. 



(20) 



Note that other forecasting papers typically assume 7r) 
of 0.3 or 0.4. 



For the CMB experiments we consider Planck [lij and 
a high-resolution version of CMBpol js^. Their spec- 
ifications are given in Table II. We assume that other 
frequency channels of Planck and CMBpol (not shown in 
the table) will clean out non-CMB sources of radiation 
perfectly. Detailed studies have shown foreground degra- 
dation of the results expected from Planck to be mild 

iii mm. At i > 3000 emission from dusty galaxies 

will be a significant source of contamination. The effect is 
expected to be more severe for temperature maps. Hence 
we restrict temperature data to I < 2000 and polarization 
data to I < 2500. 



V. ERROR FORECASTING METHOD 

The shear two-point functions depend on not just the 
dark energy parameters, but on the entire matter content 
and the primordial power spectrum. We do not assume 
these quantities to be known, but instead assume that 
CMB data are available to constrain them. The CMB 
power spectra we include in our analyses are the (un- 
lensed) Cf T , Cf E and Cf E . We do not use the lensed 
power spectra to avoid the complication of the correla- 
tion in their errors between different I values and with 
the error in the CMB-derived shear power spectrum. Us- 
ing the lensed spectra and neglecting these correlations 
can lead to overly optimistic forecasts 0. 

To calculate the expected parameter errors we make 
a first order Taylor expansion of the parameter depen- 
dence of all the CMB and cosmic shear two-point func- 
tions. In this 'linear response' approximation, given the 
expected experimental errors on the power spectra, we 
can easily calculate the expected parameter error co- 
variance matrix as the inverse of the Fisher matrix |47| . 
The linear response approximation is improved and sus- 
ceptibility to numerical error is reduced with a careful 
choice of the parameters used to span a given model space 
ElliSllllHj. As in [H] we take our (non-w(z)) set to 
be V — {u m ,u}b,L0 v ,9 s ,ZTi,k 3 Pl(kf),n s ,n' s ,y}i e }, with 
the assumption a flat universe. The first three of these 
are the densities today (in units of 1.88 x 10 _29 g/cm 3 ) 
of cold dark matter plus baryons, baryons and massive 
neutrinos. We assume two massless species and one mas- 
sive species. The next is the angular size subtended by 
the sound horizon on the last-scattering surface. The 
Thompson scattering optical depth for CMB photons, r, 
is parameterized by the redshift of reionzation z v i. The 
primordial potential power spectrum is assumed to be 
k 3 Pl(k) = kjPl,(kf)(k/k f ) n °- 1+n 's ln(k / k ^ with k f = 
0.05 Mpc -1 . The fraction of baryonic mass in Helium 
(which affects the number density of electrons) is yjjc- We 
Taylor expand about V = {0.146,0.021,0,0.6,6.3,6.4 x 
10" 11 , 1,0, 0.24}. The Hubble constant for this model is 
h = 0.655 where Ho — 100/ikmsec -1 Mpc -1 . 

We follow [H^ to calculate the errors expected in Cf T , 
Cj E and Cf E given Table 1. The errors on the unlensed 
spectra in the regime where lensing is important (deep in 
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the damping tail) are certainly underestimated because 
reconstruction of the unlensed map from the lensed map 
will add to the errors. However, this is not worrisome 
since we limit all the unlensed spectra to I < 2000. 

The contribution to the Fisher matrix from the shear- 
shear correlations is given by (e.g., [Hflj ) 

l,Bl,B2,B3,B4 

where the subscript , p denotes differentiation with re- 
spect to parameter a p . We use W to denote thet inverse 
of the total covariance matrix: 

W^S + N)- 1 . (22) 

Note that W is an easily invertible block diagonal matrix: 
yVimB,i'm'B' = W^ 5 S mm 'S B B'- The total Fisher matrix 
is given by summing F 1J with the Fisher matrix for the 
unlensed CMB data, F CMB . 

VI. RESULTS 
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Above we have described the non-w(z) parameters of 
our fiducial model, about which we consider small fluctu- 
ations. For w(z), we consider two different fiducial cases. 
Both are consant w models, one with w^d = — 1 and the 
other with wad — —0.8. About this fiducial model we 
vary all the cosmological parameters and either wo and 
w a or the five binned w^s. 

A. Dependence on Sky Coverage and Shape Noise 

First we explore the dependence of constraints on wq 
and w a as a function of sky coverage and shear weight-per 
solid angle, n/^ ms . These are both quantities over which 
the experimenter has some control. The shape of dn/dz 
is taken to be that assumed for LSST. The left panels of 
Fig. 2 show the errors on wq and w a assuming Z max = 
fOOO. The contours drop rapidly at low n/7rms but by 
^/7rms = 1000 the improvement slows dramatically. At 
such high shear weights, the dominant contribution to 
the parameter errors comes from sample variance rather 
than shear noise. 

Note that we are only varying the / s k y for the cosmic 
shear observations. We keep the / s k y for CMB observa- 
tions fixed at 1. Therefore the parameter errors reduce 

1/2 

more slowly than / sky . The improvement with increasing 
sky coverage is slowest for w a which is more dependent 
on the constraints provided by the CMB observations. 

If the non-linearities can be controlled out to I = 2000 
the errors on wq and w a improve by about 50%. Because 
the shear angular power spectrum drops with increasing 
I, the relative importance of shear noise variance to shear 
sample variance increases with increasing I. Therefore 
the contours in the right panels of Fig. 2 begin to flatten 
out at higher values of n to t / 7r ms than for the i max = 1000 
case. 



FIG. 2: Contours of constant error in wo (top panels) and w a 
(bottom panels) for cosmic shear observations of / s k y of the 
sky and shear weight per solid angle fitot/7rms' The left panels 
are for Z max = 1000 and the right panels are for Z ma x = 2000. 
These parameter error forecasts also include constraints from 
Planck. 



B. Dependence on Parameter Space and CMB 
Experiment 

Next we turn to forecasted constraints for wq and w a 
for the four cosmic shear observations combined with the 
two CMB observations. First note from Table III that 
the CMB experiments alone are incapable of constrain- 
ing these parameters. The primary CMB signal offers no 
constraint due to the angular-diameter distance degener- 
acy and CMB lensing only provides a very weak handle 
since dark energy is only important at very low redshifts, 
well below the peak of the CMB lensing kernel. The situ- 
ation improves dramatically with the introduction of the 
cosmic shear-shear correlations. 

Unlike luminosity distances to type-la supernovae, cos- 
mic shear power spectra depend on a lot of parameters 
other than Cl m , wq and w a . In fact, they depend on all 
the parameters in our full set except for the optical depth 
t. Unfortunately, the spectra are nearly featureless and 
therefore incapable of simultaneously fitting this large 
number of parameters. Thus the CMB measurements 
are key to breaking the degeneracies. 

The degeneracy-breaking provided by the CMB is 
highly effective. Forecasts for errors on the w{z) param- 
eters, assuming everything else is fixed by Planck, are 
fairly accurate as long as they include one more param- 
eter to be determined by cosmic shear data. As shown 
at the bottom of Table III, this parameter could be 8 S or 
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TABLE III: 



CMB 


Cosmic Shear 


Wo 


W a 


m v 


In fi m /i 2 


lnP* 


ns 


n's 


T 




Y P 


9, 


Planck 


none 


2.2 


3.3 


0.18 


0.0084 


0.017 


0.0078 


0.0035 


0.010 


0.0086 


0.013 


0.00016 






0.012 


0.020 




















none 


G4tt 


0.051 


0.066 


— 


0.00089 




















0.051 


0.069 


















0.00027 




S300 


0.41 


0.49 


0.12 


0.0068 


0.016 


0.0066 


0.0034 


0.0091 


0.0076 


0.011 


0.00014 


Planck 


S1000 


0.28 


0.34 


0.090 


0.0059 


0.015 


0.0060 


0.0033 


0.0090 


0.0069 


0.011 


0.00014 




G2tt 


0.076 


0.11 


0.052 


0.0039 


0.014 


0.0046 


0.0029 


0.0083 


0.0059 


0.0094 


0.00012 




G4n 


0.056 


0.087 


0.045 


0.0035 


0.013 


0.0044 


0.0025 


0.0077 


0.0058 


0.0087 


0.00012 






0.049 


0.063 




0.00083 




















0.021 


0.031 


















9.4 x 10~ 5 






0.050 


0.064 




0.00086 














9.7 x 10 -5 


Planck 


G4tt 


0.050 


0.065 


0.018 


0.0015 














9.8 x 10~ 5 






0.050 


0.069 


0.018 


0.0028 


0.0093 




— 


0.0055 


— 


— 


0.00011 






U.Uol 


U.U ( U 


U.UzU 


U.UUoo 


n nil 


U.UUzU 




U.UU f U 






U.UUU11 






0.052 


0.074 


0.023 


0.0033 


0.011 


0.0020 


0.0019 


0.0070 






0.00011 




none 


0.79 


1.1 


0.048 


0.0040 


0.017 


0.0031 


0.0018 


0.0097 


0.0028 


0.0048 


4.8 x 10" B 




S300 


0.18 


0.24 


0.047 


0.0038 


0.012 


0.0029 


0.0018 


0.0073 


0.0028 


0.0047 


4.6 x 10~ 5 


CMBpol 


S1000 


0.13 


0.17 


0.045 


0.0036 


0.011 


0.0028 


0.0018 


0.0069 


0.0028 


0.0046 


4.6 x 10" 5 




G2tt 


0.064 


0.089 


0.031 


0.0026 


0.0099 


0.0026 


0.0016 


0.0060 


0.0027 


0.0044 


4.2 x 10" 5 




G4tt 


0.049 


0.070 


0.027 


0.0023 


0.0093 


0.0025 


0.0015 


0.0057 


0.0026 


0.0043 


4.1 x 10" 5 



NOTES. — Forecasts of 1-<t errors for the combination of tomographic shear-shear two-point 
functions with Planck and CMBpol. A ' — ' indicates a parameter held fixed. 



Q m h 2 . However, assuming that Planck determines both 
9 S and fl m h 2 perfectly well results in overly optimistic 
forecasts for wq and w a . This is a bit surprising since the 
CMB does such an exquisite job on both 9 S and Vl m h 2 . 
The reason for it is that the shear power spectra are 
highly sensitive to the shape parameter, U, m h which is 
a combination of 6 S and fl m h 2 . Even very small uncer- 
tainties in Q m h can significantly increase the errors on 
Wq and w a - 

Replacing Planck with four-year WMAP, the error 
forecasts for wq and w a from G4-7T weaken from 0.056 
to 0.064 and from 0.087 to 0.11 respectively. Replacing 
Planck with CMBpol, as we see in Table III, strengthens 
the constraint to the limit of all parameters (other than 
wo, w a and either 9 S or fl m h 2 ) being fixed by the CMB. 



C. Neutrino Mass 

Replacing a massless neutrino species with a massive 
one increases the energy density and therefore the ex- 
pansion rate, suppressing the growth of structure. The 
suppression of the power spectrum is scale dependent and 
the relevant length scale is the Jeans length for neutrinos 
I55L which decreases with time as the neutrino 
thermal velocity decreases. This suppression of growth 
is ameliorated at scales larger than the Jeans length at 



matter-radiation equality, where the neutrinos can clus- 
ter. Neutrinos never cluster at scales smaller than the 
Jeans length today. The net result is no effect on scales 
larger than the Jeans length at matter-radiation equal- 
ity and a scale-indepenent suppression of power on scales 
smaller than the Jeans length today. Both the total sup- 
pression of power and the Jeans length today depend on 
the neutrino mass. 

The observed galaxy power spectrum (which is propor- 
tional to the matter power spectrum at sufficiently large 
scales) , combined with CMB observations can be used to 
put constraints on m„ [5?|. At present such an analysis 
yields an upper bound on m v of ~ 0.3 eV |58| . 

Since the shear power spectra depend on the history 
of the matter power spectrum they can also be used to 
determine neutrino masses. Recently, [241 ] showed that an 
experiment with the specifications of CMBpol in Table II 
could determine the mass of a relic neutrino species via its 
effect on the lensing of the cosmic microwave background 
with an uncertainty of 0.044eV. This analysis assumed a 
time-independent w. We see in Table III that allowing 
for variation of w a weakens the neutrino constraint only 
slightly to 0.048 eV. A reasonable prior on wq and w a 
would eliminate this degradation. 

Adding constraints from shear-shear correlations to 
CMBpol only improves the neutrino mass constraints by 
less than a factor of 2. We will see more dramatic im- 
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provement with the introduction of the ratio statistic in 
the next section. Shear-shear correlations make a more 
dramatic improvement for m v when added to Planck. 
The combination of Planck and G2-7T does nearly as well 
as CMBpol alone. A similar result was found in |22u23| . 
although with the assumption of w = — 1. 



D. Primordial Power Spectrum Parameters 

One of the virtues of inflation is the robustness and 
simplicity of its predictions. On the other hand, this ro- 
bustness and simplicity make it difficult to probe obser- 
vationally in detail. One of the few observational handles 
we have is the shape of the primordial power spectrum. 
High-precision determination of ng and (even better) de- 
termination of non-zero n' s = dns/d In k would be highly 
valuable information about inflation or whatever mecha- 
nism generated the initital fluctuations. 

Recently the scalar power spectrum has been a sub- 
ject of great interest due to claims of non-zero n' s . The 
evidence from CMB data alone is very weak: n' s = 
dn s /d\nk = -0.055 ±0.038 58]. Combining CMB data 
with the Croft et al. matter power spectrum inferred 
from high-resolution observations of the Ly Q forest re- 
sults in n' s = -0.031 ± 0.017 [58|. Other authors [U 
working with the same datasets though have since found 
looser constraints when they marginalize over the mean 
ionizing flux as a function of redshift. The importance of 
marginalizing over this parameter, which leads to a large 
degeneracy between spectral index and amplitude, was 
pointed out in |6*c|. 

Progress in studying the scalar power spectrum is 
likely to come mostly from higher-resolution CMB exper- 
iments. As we can see in Table III (and in [24]]) Planck 
will improve the constraints on rig from current values 
by a factor of 10. If the error bars shrink around the 
current best-fit values, this improved precision will have 
profound implications since such large departures from 
scale-invariance are not expected in inflationary models. 

The difference ng — 1 is first order in the slow-roll pa- 
rameters and n' s is second order. Thus we expect \n' s \ 
to be on the order of (ns — l) 2 - If 1 — ng = 0.03, which 
is perfectly consistent with present data, then we would 
expect to see non-zero n' s at about the 10 -3 level. It 
would be tremendously exciting to actually confirm this 
expectation. Note though that there are models with 
| rig | ~ ns — 1 [6l| so this is not a firm test. Nevertheless, 
finding cr(n' s ) » ng — 1 would be difficult to reconcile 
with inflation |6l|. 

Confirming the expectation that ng — 1 is on the order 
of (ng — l) 2 is likely to take the precision of CMBpol, 
as pointed out in |24j . Unfortunately, even if 1 — ng is 
as large as 0.06, CMBpol alone would still only achieve 
a two-sigma detection of n' s — (ns — l) 2 . Further, the 
addition of cosmic shear to CMBpol does not improve 
the constraints on n' s by much, although it can improve 
cr(n'g) by a factor of five over what's possible with WMAP 



data alone [62| . 

It is very hard to imagine what could improve the de- 
termination of n' s further • Extending dynamic range 
further only improves the error logarthmically, and then 
only if power at the smaller scale is measured with pre- 
cision comparable to that at the larger scale. It is highly 
unlikely that systematics on any small-scale measure- 
ment could be controlled well enough to result in con- 
straints on the primordial power spectrum amplitude at 
the 10 -3 level. 

Finally, note that the amplitude of the primordial 
power spectrum is better-determined with the addi- 
tion of cosmic shear. And since the CMB determines 
P$ exp(— 2r) to very high accuracy, the r constraint im- 
proves too |3 ■ 

E. Dark Energy Tomography 

Now we turn to the results for the discrete parametriza- 
tion of w(z), shown in Fig. 3 and Fig. 4. We see that the 
errors grow with increasing redshift as expected due to 
the decreasing importance of dark energy with increasing 
redshift. For w = — 1, matter-dark energy equality is at 
z eq ~ 0.25. The percentages of the total density in dark 
energy are 25%, 7% and 3% at z = 1, z = 2 and z = 3 
respectively. For w — —0.8 the drop off is slower so the 
rise in the errors is slower. The highest-redshift bin is 
the broadest bin, stretching from z = 1.6 all the way to 
z = 1100; this width is responsible for the smaller error 
in the highest redshift bin. 

We see from Fig. 3 and Fig. 4 that weak lensing ob- 
servations of very large amounts of sky will allow for a 
reconstruction of w(z) up to z — 2, with errors in all five 
bins (for the G47T, Wfid = —0.8 case) less than 0.2. 

The smaller sky coverage of the space-based surveys 
leads to much worse constraints. Although here we must 
remind the reader that we are focusing on science that 
can be done with information on the larger angular scales. 
Space-based observations will have advantages with re- 
spect to ground-based observations at smaller angular 
scales. 

The errors in u>i are highly correlated. It is useful for 
us to find an uncorrelated linear combination of W{ when 
we fit the large scale modes at the second stage. We write 
w a = RaiWi and find R a i such that 

(Aw a Aw p) = F~p = E 2 J afS 

5 

= 2J RaiF^Rjp (23) 

where F~} is the (diagonal) Fisher matrix for the trans- 
formed parameters and the transformation matrix R sat- 
isfies R a iRfp = S a p. TABLE IV shows 1-a errors for the 
amplitudes of the eigenmodes of Wi . 

After the first stage of analysis, restricted to I > 40 and 
ignoring dark energy fluctuations, one would proceed to 
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w rid=-10 
G47r + Planck 



w fld =-0.8 
G47T + Planck 



1 1 

w rid=-10 
G47T + CMBpol 



w rid =-0.8 
G47T + CMBpol 
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Z 



0.4 0.8 1.2 1. 
Z 



FIG. 3: Error boxes indicate expected 1-a error in each Wi 
for G47r combined with either Planck (left panels) or CMBpol 
(right panels) for w^d = —1 (top panels) or -0.8 (bottom 
panels) . 







El 


E2 | E3 | E4 | E5 


0.021 


0.046 1 0.088 1 0.13 1 0.28 







TABLE IV: Forecasts for 1-a errors on the Wi eigenmode am- 
plitudes for G47T and CMBpol assuming wad = — 1. 



a second stage which would extend to larger scales and 
take the fluctuations into account. The second stage of 
analysis would use a smoothed version of the eigenmodes 
derived from the first stage of analysis as the basis for 
the parametrization of w(z). We can test the consistency 
from what we found at the first stage, though we do not 
expect noticeable improvement at the second stage. Be- 
cause the large-scale and small-scale data are uncorre- 
cted, the small-scale data could simply be included as a 
prior determined in the first stage of the analysis. Iclud- 
ing the small-scale data in this way may be of practical 
use, in terms of making the data analysis tractable, if 
the parameter space needs to be broadened to include 
the sound-speed of the dark energy. For a study of how 
well the sound speed can be constrained see 0. 



VII. THE RATIO STATISTIC 



, — ,w 



.w M =-1.0 
G4tr+Planck 



■1.0 



SI 000 + Plane* 



1 1 

, Wfid =-1.0 

G27T+Planck 



w rid=-l-3 
S300+P anck 



0.4 0.8 1.2 1.6 0.4 0.8 1.2 1.6 

Z Z 

FIG. 4: Error boxes indicate expected 1-a error in each Wi 
for the combination of experiments specified in each panel. 



We have restricted our use of the shear maps to I < 
1000 to reduce our sensitivity to theoretical errors in the 
calculation of P(k, z). The density power spectrum is 
difficult to calculate on small scales where the influence 
of baryons becomes important. 

However, discarding the data at I > 1000 is highly 
unfortunate since there is tremendous signal-to-noise on 
these angular scales, as can be seen in Fig. 1. Recently 
Jain & Taylor |2(j proposed a ratio statistic that uses 
this small-scale data to constrain cosmology, but that is 
not biased by errors in the small-scale power spectrum. 
Here we combine the ratio statistic with the shear-shear 
power spectra at I < 1000 and find it provides highly 
complementary constraints on the bins of w. 

Further details on how the Jain & Taylor idea l20l 
might be implemented are given in Bernstein & Jain [21| • 
They write the observed eccentricity as 



..obs 



I 



di{yik)gik + e fc . 



(24) 



The distortion field di is the distortion field for sources 
at infinity due to mass inside redshift shell I. The gik is 
a purely geometric term that corrects the distortion for 
a source at finite distance r^. If the distortion field is 
dominated by a single lensing object then we can write 



^obs 



= du(x k )gi*k + / y di(x. k )gik + e fc (25) 



ijti* 
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I I I I I I I I I 



J-1.4 -1.2 -1 -0.8 -0.6 J-1.4 -1.2 -1 



-0.8 -0.6 




FIG. 5: Forecasts of l-cr error contours for pairs of Wi assum- 
ing G4-7T and CMBpol. The thin solid contour is for shear- 
shear correlations alone, the dashed contour is for the ratio 
statistic alone and the thick contour is for their combination. 



which is just rewriting the above, but with the under- 
standing that the term with the sum is now small com- 
pared to the others. With the assumption of a shape 
for the distortion field, di*(x.k), one can determine guk 
up to a normalizing constant, as well as gu k ' up to the 
same normalizing constant. Thus Jain and Taylor pro- 
posed using the ratio of such determinations to give a 
purely geometric quantity. This is their ratio statistic. 
Its expectaion value is 



-Rfefc' — Qik/gik'- 



(26) 



A rough idea of how well R kk , can be determined is 
given by considering the typical amplitude of the shear 
signal and the area over which it can be measured, the 
latter affecting how many galaxies can be used to beat 
down shape noise. This leads to 



^kk',jj' 



AR kk' 



Ml 



kk' 



Sk 



4') 



5k'r 



(27) 



where j k is the average amplitude of the tangential shear 
of galaxies in source bin k due to lenses in lens slice /. 
The noise variance (n k ) is given by 



7 2 

/rms 

JfkJi 



(28) 



where /; denotes the fraction of the sky covered by tem- 
plates for lenses in lens redshift slice I. 

To proceed further we need estimates for the shear 
template amplitudes -y k and the fraction of the sky in 
the observed regions covered by the templates, Jain & 
Taylor 20] point out that the number density of clusters 
with mass in a logarithmic intermal about m = 10 14 M 
out to z = 1 is 0.01/arcmin 2 . Further, for a galaxy clus- 
ter of this mass at z = 0.2 with lens source at z = 1 the 
tangential shear signal is nearly constant out to 1' with 
an area-weighted average amplitude of about 0.05 [&^| . 
We therefore set j k — 0.05 and define the area useable 
for estimating the template amplitude to be that within 
the 1' radius. With this template area per object and 
considering objects with masses in the decade surround- 
ing m = 10 14 Mq we get a fraction of the observing area 
covered with templates of about 0.06. Dividing this into 
three lens slices (centered on z = 0.2, z = 0.6 and z — 1) 
we get ji = 0.02. 

Jain and Tavlorpol] used a full halo model calculation 
of the signals expected and area that could be used for 
measuring them in each lens slice for each mass of clus- 
ter in order to forecast errors on their ratio statistic, and 
therefore on dark energy parameters. They also used 
a version of the simple forecast above that was able to 
reproduce the results of their more detailed calculation. 
Our simple forecast for errors on the ratio statistic is 
much more conservative than theirs. The main difference 
is that they assumed an area with a radius of 3' about 
each cluster could be used for measuring the shear tem- 
plate amplitude. We chose the smaller radius to reduce 
contamination from projection. 

Our results including lcnsing cosmography are more 
speculative than our results from just CMB and the shear 
two-point functions. Projection effects may be impor- 
tant, even with our restriction to 1'. Although [21| in- 
cluded projection effects, they did so with a halo model 
calculation which may be inadequate for this purpose 
given its assumption of spherical halos. Further, the con- 
struction of templates from the gala xy distribution may 
not be as robust as assumed in |2l| . Numerical simula- 
tions could address the projection problem [l^. Finally, 
unlike the sample- variance dominated shear-shear corre- 
lation statistics, the errors we forecast on the ratio statis- 
tic depend strongly on the achievement of the low level 
of shape noise that we assume. 

The Fisher matrix for the ratio statistic is given by 



F, 



i>i>' 



EE 

/ kk' 



R 



R 



kk',p „l _i 

^kk'.kk' 



Ri 



kk' ,p' 



kk' 



R 



i 

kk' 



(29) 



where we have approximated the C matrix to be inde- 
pendent of cosmological parameters. We calculate this 
Fisher matrix summing over the three lens planes and 
all eight cosmic shear source planes. Because the ratio 
statistic is influenced by data on smaller angular scales 
than the shear-shear statistics (since there we restricted 
to I < 1000) the errors are independent and the Fisher 
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TABLE V: 



CMB 


1 nctnic s nociT' 


Wo 


W a 




In O h 2 
ill i &rn i v 


In P<$ 


ns 


n s 




In a; 


i p 


Q 

"s 




none 


2.2 


3.3 


0.18 


0.0084 


0.017 


0.0078 


0.0035 


0.010 


0.0086 


0.013 


0.00016 




S300 


0.16 


0.26 


0.074 


0.0061 


0.015 


0.0065 


0.0033 


0.0090 


0.0073 


0.011 


0.00014 


Planck 


S1000 


0.10 


0.16 


0.051 


0.0048 


0.015 


0.0058 


0.0031 


0.0089 


0.0067 


0.011 


0.00014 
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n m q 


n no^n 

u.uuoo 


o oi a 


n 00^7 

U.UUO 1 


n 0090 


0080 


00^8 
u.uuoo 


00fi7 

U.UUU 1 


0001 9 




G4tt 


0.022 


0.040 


0.014 


0.0030 


0.013 


0.0033 


0.0017 


0.0074 


0.0056 


0.0059 


0.00011 




none 


0.79 


1.1 


0.048 


0.0040 


0.017 


0.0031 


0.0018 


0.0097 


0.0028 


0.0048 


4.8 x 10 -5 




S300 


0.12 


0.20 


0.044 


0.0036 


0.012 


0.0029 


0.0018 


0.0073 


0.0028 


0.0046 


4.6 x 10~ 5 


CMBpol 


S1000 


0.082 


0.13 


0.037 


0.0033 


0.011 


0.0028 


0.0018 


0.0066 


0.0028 


0.0046 


4.5 x 10 -5 




G2tt 


0.028 


0.047 


0.017 


0.0022 


0.0086 


0.0024 


0.0014 


0.0053 


0.0026 


0.0037 


3.9 x 10 -5 




G4tt 


0.021 


0.035 


0.013 


0.0021 


0.0080 


0.0021 


0.0013 


0.0049 


0.0025 


0.0033 


3.6 x 10" 5 



NOTES. — Forecasts of 1-a errors for the combination of tomographic shear-shear two-point 
functions and the ratio statistic with Planck (upper table) and CMBpol (lower table). 
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FIG. 6: Forecasts of 1-a errors in Wi for the i = 1 to 5 redshift 
bins for CMBpol combined with the all-sky weak lensing ex- 
periment G47T. The thin solid error boxes are for shear-shear 
correlations alone, the dash error boxes are for the ratio statis- 
tic alone and the thick solid boxes are for the combination. 



matrices simply add: 

Fpp' — Fppi + Fp^pi + Fp p , . (30) 

Despite our more conservative treatment, we still find 
the ratio statistic to provide a powerful constraint on 
dark energy especially when combined with the shear- 
shear correlations investigated above. Parameter error 
forecasts for the two-parameter dark energy model from 



combining shear-shear correlations and the ratio statis- 
tic with CMB data are shown in Table V. We see that 
the addition of the ratio statistic reduces errors in the 
dark energy parameters by about a factor of two. Con- 
straints on neutrino mass also improve by factors of two 
and three. 

Even more dramatic improvements occur in the Wi 
space. The combination of shear-shear correlations and 
the ratio statistic dramatically reduces the degeneracies 
of these parameters. Fig. 5 shows that with shear-shear 
correlations the errors in Wi and Wi + i are anti-correlated 
while with the ratio statistic the errors are positively cor- 
related. The result for the combination is a large reduc- 
tion of the error for each Wi as can be seen in Fig. 6. 
As is shown in Fig. 6, G4-7T plus CMBpol can reconstruct 
w(z) with uncertainty a(wi) < 0.1 up to z ~ 2. 

The anti-correlation from shear-shear correlations is 
expected since a positive variation of Wi can be com- 
pensated for by a negative variation of iWj+x. But the 
complicated ratio statistic leads to the unusual postive 
correlation, where a positive variation of Wi can be par- 
tially compensated for by a positive variation in tUf+i. 



VIII. DISCUSSION AND CONCLUSIONS 

The statistical properties of the large-scale density field 
are sensitive to the dark energy equation of state. Tomo- 
graphic cosmic shear observations on large angular scales 
offer a theoretically robust means of taking advantage of 
this dependence. While the statistical properties are also 
sensitive to other matter components and the primordial 
power spectra, in combination with CMB observations 
dark energy parameters can still be precisely measured 
despite these other degrees of freedom. 

We have seen that very tight constraints can be 
achieved on the parameters of an 11-dimensional space, 
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including the equation-of-state parameter and its time 
variation, the mass of the relic neutrinos, and the shape 
and amplitude of the primordial power spectrum. Cosmic 
shear surveys are critical to the dark energy constraints, 
and constributc significantly to the other constraints, es- 
pecially those on neutrino mass. 

To explore the constraints possible on the history of the 
dark energy equation-of-state parameter we introduced a 
discrete parametrization oiw(z). The combination of all- 
sky weak lensing observations on large angular scales and 
all-sky CMB observations can constrain w in each of five 
bands (four evenly spaced from to 1.6 and one from 1.6 
to last-scattering) to better than 0.2. 

Such precision on w out to redshifts where the dark 
energy density is subdominant may prove to be highly 
illuminating as to the nature of the dark energy. Some 
attractor solutions predict a sharp transition of w from 
m > to m = - lata redshift where the dark energy 
density is still subdominant |65l I66L |6j| ■ If such a transi- 
tion happens at z < 2, then we can probe that transition 
point with the binned «;,. 

While we have focused on what is possible with the 
statistical properties that can be most robustly predicted 
from theory, there is much signal to noise at smaller an- 
gular scales. Further work is needed to make full use of 
these data where non-linear corrections are large. We 
studied one means of including information from small- 
scales, the ratio statistic of [2(J], and found it resulted 
in large reductions in forecasted parameter errors. The 
m v constraint for G47r+Planck improved by a factor of 
three to the very interesting level of 0.014 eV. The G4-7T 



+ CMBpol constraints on the Wi improved even more 
dramatically in the three lowest redshift bins. 

Many challenges remain to the achievement of the pa- 
rameter errors we forecast here. We have assumed that 
residuals from fitting anisotropic point-spread functions, 
distortions introduced by the camera, photometric red- 
shift errors, redshift-dependent calibration errors |68| and 
intrinsic alignments |69j, l70l l7ll 1 721 l73l |74j are negligible. 
The validity of these assumptions depends on the skill 
of the experimenters and analysts and to some degree on 
the kindness of nature. Study of the impact of these non- 
idealities and how to tackle them is on-going. Calibration 
errors of 5%, achievable with current data analysis meth- 
ods [6^ , lead to only small changes in forecasted parame- 
ter e~rs M- Intrinsic alignment contamination can be 
avoided by using photometrically-determined redshifts to 
exclude spatially-close galaxy pairs[8^ or, more radically, 
by restricting analyis to shear cross-spectra and ignoring 
the auto spectra |75j . Avoiding significant biases in pho- 
tometric redshifts will require calibration with larg e, but 
achievable, numbers of spectroscopic redshifts |2lj . 
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